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Abstract 

The angular distributions of the baryon-antibaryon low-mass enhancements seen in the charmless 
three-body baryonic B decays B + — > ppK + , B° — > ppKg, and B° — > pAir~ are reported. Searches 
for the pentaquarks + and B ++ in the relevant decay modes and possible glueball states in the 
pp systems are presented. The analysis is based on a 140 fb _1 data sample recorded on the T(4S) 
resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. 

PACS numbers: 13.25.Hw, 13.60.Rj 
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Observations of several baryonic three-body B decays have been reported recently |l| 



mass spectra toward threshold, as originally conjectured in Refs. |4|, |5( and elaborated more 
recently in Refs. jfj 0, H[. The same peaking behavior near threshold has been found in 
J I ip decays |9| as well, indicating that this might be a universal phenomenon. The possible 
explanations include the presence of intermediate gluonic states or side effects of the quark 
fragmentation process. Alternatively, the dynamical picture can be replaced by an effective 
range analysis with a baryon form factor. To distinguish among the above production 
mechanism hypotheses, we study the threshold enhancements by examining the angular 
distributions in the helicity frame for the ppK + , ppK s and pAir~ ^(J modes. Also, we 
update the mass spectra from our previous studies. 

We use a 140 fb _1 data sample, consisting of 152 xlO 6 BB pairs, collected by the Belle 
detector at the KEKB asymmetric energy e + e~ (3.5 on 8 GeV) collider fllj. The Belle 
detector is a large solid angle magnetic spectrometer that consists of a three layer silicon 
vertex detector (SVD), a 50 layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time of flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return 
located outside of the coil is instrumented to detect K\ mesons and to identify muons. The 
detector is described in detail elsewhere |l2j |. 

The event selection criteria are based on the information obtained from the tracking sys- 
tem (SVD+CDC) and the hadron identification system (CDC+ACC+TOF). All primary 
charged tracks are required to satisfy track quality criteria based on the track impact param- 
eters relative to the interaction point (IP). The deviations from the IP position are required 
to be within ±1 cm in the transverse (x-y) plane, and within ±3 cm in the z direction, 
where the z axis is opposite the positron beam line. For each track, the likelihood values 
L p , Lk-i and L n that it is a proton, kaon, or pion, respectively, are determined from the 
information provided by the hadron identification system. The track is identified as a pro- 
ton if L P /(L P + Lk) > 0.6 and L P /(L P + L n ) > 0.6, or as a kaon if Lk/(Lk + L n ) > 0.6, 
or as a pion if L n /(Lx + L n ) > 0.6. Candidate K s mesons are reconstructed from pairs of 
oppositely charged tracks (both treated as pions) having an invariant mass consistent with 
the K s nominal mass, as well as a displaced vertex and flight direction consistent with an 
origin at the IP. A candidate A baryon is reconstructed from a pair of oppositely charged 
tracks — treated as a proton and negative pion — whose invariant mass is consistent with the 
nominal A baryon mass. The proton-like daughter is required to satisfy L p / (L p + L n ) > 0.6. 

Candidate B mesons are reconstructed from the related final state particles for the B + — > 
ppK + , B° — > ppK s , and B° — > pAir~ modes. We use two kinematic variables in the center 
of mass (CM) frame to identify the reconstructed B meson candidates: the beam energy 



constrained mass Mb c = y E^ eam — p B , and the energy difference AE = E B — E hcSLm , where 
-E-beam is the beam energy, and ps and E B are the momentum and energy, respectively, of 
the reconstructed B meson. The fit region is defined as 5.20 GeV/c 2 < Mb c < 5.29 GeV/c 2 
and -0.1 GeV < AE < 0.2 GeV. From a GEANT based Monte Carlo (MC) simulation, the 
signal peaks in the subregion 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and \ AE\ < 0.05 GeV. The 
lower bound of AE is chosen to exclude possible contamination from so-called "cross-feed" 
baryonic B decays. 

The background in the fit region solely arises from the continuum e + e~ — > qq (q = 
u, d, s, c) process. We suppress the jet-like continuum background events relative to the 




4 



more spherical BB signal events using a Fisher discriminant [13[ that combines seven event 
shape variables, as described in Ref. [14|. Probability density functions (PDFs) for the 
Fisher discriminant and the cosine of the angle between the B flight direction and the beam 
direction in the T(4S) rest frame are combined to form the signal (background) likelihood C s 
{Cb). The signal PDFs are determined using signal MC simulation; the background PDFs 
are obtained from the side-band data with M bc < 5.26 GeV/c 2 . We require the likelihood 
ratio 1Z = C s / (C s + £&) to be greater than 0.7, 0.75, and 0.8 for the ppK + , ppK$, and pAn~ 
modes, respectively. These selection criteria are determined by optimization of n a j \Jn s + n b , 
where n s and n& denote the expected numbers of signal and background events, respectively. 
We use the branching fractions from our previous measurements 0, 0, HI in the calculation 
of n s . If there are multiple B candidates in one event, we select the one with the best x 2 
value from the B meson vertex fit, in which only the primary charged tracks are used. Based 
on previous studies [HQ, 0], we require the invariant mass of the baryon pair to be less than 
2.85 GeV/c 2 for the study of the threshold enhancement effect that follows. 

The Mbc distributions (with \AE\ < 0.05 GeV), and the AE distributions (with Mb c > 
5.27 GeV/c 2 ) for the ppK + , ppKg and pAir~ modes are shown in Fig.^ We use an unbinned 
likelihood fit that maximize the likelihood function, 

p -(N s +N b ) N 

L = I] [N s P s (M hCi7 AE,) + N b P b (M hCi , AEj] , 

to estimate the signal yield; here P s (Pb) denotes the signal (background) PDF, N is the 
number of events in the fit, and N s and N b are free parameters representing the number of 
signal and background events, respectively. 

For the signal PDF, we use the product of a Gaussian in M bc and a double Gaussian 
in AE. We fix the parameters of these functions to values determined by MC simulation 
[lfij ]. The continuum background PDF is taken as the product of shapes in M bc and AE, 
which are assumed to be uncorrelated. These shapes are obtained from sideband events, 
with 0.1 GeV < AE < 0.2 GeV for the M bc function and with 5.20 GeV/c 2 < M bc < 
5.26 GeV/c 2 for the AE function, and are confirmed with a continuum MC sample. We 
use the following parametrization first used by the ARGUS collaboration 16], /(M bc ) oc 
M bc \/1 — x 2 exp[— £(1 — x 2 )], to model the M bc background, with x given by M hc /E heam and 
£ as a fit parameter. The AE background shape is modeled by a first order polynomial 
whose slope is a fit parameter. The projections of the fit results are shown in Fig. by solid 
curves. The fit yields are 217 ± 17, 28.6 tjjjj, and 48.8 t?j for the ppK + , ppK° s , and pAn~ 
modes, respectively. 

We study the proton angular distribution of the baryon-antibaryon pair system in its 
helicity frame. The angle 8 P is defined as the angle between the proton direction and the 
opposite of the light meson direction in the baryon-antibaryon pair rest frame. Note that 
after charge conjugation, the angle is determined by p and K + (or p and K~) for the ppK + 
mode. Fig. |2fa)-(c) shows the branching fractions as a function of cos8 p . The error bars 
include the statistical uncertainty from the fit and the systematic uncertainty. It is clear 
that the fragmentation process is favored for the ppK + mode. Protons are emitted along 
the K~ direction most of the time, which can be explained by a parent b — > s penguin 
transition followed by su fragmentation into the final state. The cos# p distribution of the 
ppK^ mode does not have enough statistics to support or refute this interpretation, although 
it seems to be peaked towards ±1 since the flavor information is not applied in this case. 
The distribution for the pA7r~ mode is quite flat, in contrast to that of the ppK + mode, 
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although both presumably share a common origin in the b — > s transition. In fact, this 
parentage suggests that it would be more suitable to draw the proton angular distribution 
of the pir~ pair relative to the A direction; this is shown in Fig. HJcl). It is evident that 
the above interpretation is supported: the proton tends to emerge parallel to the A baryon. 
As a cross check, the distribution of cos8 p for background events in the ppK + sample is 
shown in Fig.^e). (Similar distributions are obtained for the backgrounds of the ppK° s and 
pKir~ modes.) The background has a 1 + cos 2 9 P distribution, which can be explained as 
arising from the random combination of two high momentum particles from the qq jets. The 
fragmentation signature is not seen in the B + — > J/ipK + mode, where the J/ip meson decays 
to app pair. For J/ip candidates with invariant mass in the range 3.07 GeV/c 2 < M pp < 3.11 
GeV/c 2 , the distribution of cos6 p is flat, as shown in Fig. 

The differential branching fraction as a function of the baryon pair invariant mass is 
shown in Fig. |3| Here, the efficiency as a function of baryon pair mass for each signal mode 
is determined by MC simulation, with the events distributed uniformly in phase space. The 
regions 2.850 GeV/c 2 < M pP < 3.128 GeV/c 2 and 3.315 GeV/c 2 < M pP < 3.735 GeV/c 2 are 
excluded to remove background from B decay modes containing an r) c , J/ip, ip', Xco, or Xci 
meson. The width of the low mass enhancement in each distribution of Fig. |3] depends on 
the signal mode. A different narrow width is seen also in the newly discovered B + — > AAK + 
decay 

Systematic uncertainties are determined using high statistics control data samples. For 
proton identification, we use a A ^ pn~ sample, while for K/tt identification we use a 
D* + — > D°7i + , D° — > K~7i + sample. Tracking efficiency is measured with fully and partially 
reconstructed D* samples. The Kg reconstruction efficiency is determined from a D~ — > 
K s tt~ sample. The A and K° s reconstruction efficiencies have the same uncertainty due to 
off-IP tracks if the uncertainty of the daughter proton identification criterion is not taken 
into account. The 1Z continuum suppression uncertainty is estimated from b — > c control 
samples with similar final states. Based on these studies, we assign a 1% error for each 
track, 3% for each proton identification, 2% for each kaon/pion identification, 5% for Kg 
and A off-IP reconstruction and 6% for the 1Z selection. 

A systematic uncertainty of 4% in the fit yield is determined by varying the parameters of 
the signal and background PDFs. The MC statistical uncertainty and binning of the baryon 
pair mass contribute a 2% error in the branching fraction determination. The error on the 
number of BB pairs is 1%, where the assumption is made that the branching fractions of 
T(4S) to neutral and charged BB pairs are equal. 

We first sum the correlated errors linearly and then combine with the uncorrelated ones 
in quadrature. The total systematic uncertainties are 11%, 12%, and 12% for the ppK + , 
ppK s , and pAir~ modes, respectively. 

The newly observed narrow pentaquark state, + |l8[, can decay into pK$. We perform 
a search in our data sample by requiring 1.53 GeV/c 2 < M pK o < 1.55 GeV/c 2 . The M bc 
and AE projection plots in Fig. 0] show no evidence for a pentaquark signal. Since there 
are few events in the fit window, we fix the background shapes from side-band data. We 
use the fit results to estimate the expected background and compare this with the observed 
number of events in the signal region to set the upper limit on the yield 0, The 
systematic uncertainty is included in this estimation. The upper limit yield is determined to 
be 3.9 at the 90% confidence level. The related upper limit product of branching fractions is 
B(B° -> Q + p) x B(Q + -> pK%) < 2.3 x 10" 7 at the 90% confidence level. We also perform a 
search for ++ , which can decay to pK + . Because there are only theoretical conjectures for 
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the existence of such a state, we examine the wider mass region of 1.6 GeV/c 2 < M p x+ < 1.8 
GeV/c 2 . We find no evidence for signal. Assuming this state is narrow and centered near 
1.71 GeV/c 2 , the upper limit yield is 3.3 events at the 90% confidence level. The related 
upper limit product of branching fractions is B(B + — > Q ++ p) x i3(0 ++ — > pK + ) < 9.1 x 10 -8 
at the 90% confidence level. 

One theoretical conjecture suggests that a glueball resonance with mass near 2.3 
GeV/c 2 may explain the M pp threshold peaking behavior for the ppK + mode. Since the M pp 
mass resolution is about 10 MeV/c 2 , we scan through the 2.2 GeV/c 2 < M pp < 2.4 GeV/c 2 
mass region with a 20 MeV/c 2 wide window. The highest upper limit yield is found to be 
18.9. We use this data set to set an upper limit on the product of branching fractions of 
B(B + — > glueball K + ) x £> (glueball — > pp) < 4.1 x 10~ 7 at the 90% confidence level for 
a possible narrow glueball state with mass in the 2.2 - 2.4 GeV/c 2 range. The theoretical 
prediction is around 1 x 10 -6 . 

In summary, using 152 xl0 6 BB events, we measure the angular and invariant mass 
distributions of the baryon-antibaryon pair system near threshold for the ppK + , ppK® and 
pA7r~ baryonic B decay modes. The quark fragmentation process is supported, but the 
gluonic picture is disfavored. Searches for a B meson decaying into pentaquark B + or a 
glueball in the above related modes give null results. We set stringent upper limits on the 
product of the decay branching fractions. 
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FIG. 1: Distributions of AE and M^ c , respectively, for (a) and (b) ppK + , (c) and (d) ppK$, and 
(e) and (f) pA.7r~ modes with baryon-antibaryon pair mass less than 2.85 GeV/c 2 . The blue, red 
and dashed lines represent the combined fit result, fitted signal and fitted background, respectively. 
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FIG. 2: Branching fraction vs. proton helicity angle in the baryon-antibaryon pair system for (a) 
ppK + , (b) ppKg, and (c) pAir~ modes, (d) The proton angular distribution of the pir~ system 
against the A direction in the pAtt~ mode, (e) Background yield of ppK + in the fit. (f ) Comparison 
with the J/ip mass region. 
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FIG. 3: Differential branching fraction for (a) ppK + , (b) ppK^, and (c) pA7r~ modes as a function 
of baryon-antibaryon pair mass. The shaded distribution shows the expectation from a phase-space 
MC simulation with area scaled to the signal yield. A charmonium veto has been applied in (a) 
and (b). 
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FIG. 4: Distributions of M bc and A£ for the ppK° s mode with 1.53 GeV/c 2 < M pi ^o < 1.55 
GeV/c 2 . The curves represent the fit projections. 
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